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COUPLED PARAXIAL WAVE EQUATIONS IN RANDOM MEDIA 
IN THE WHITE-NOISE REGIME^ 

By Josselin Garnier and Knut S0lna^ 
Universite Paris VII and University of California at Irvine 

In this paper the reflection and transmission of waves by a three- 
dimensional random medium are studied in a white-noise and parax- 
ial regime. The limit system derives from the acoustic wave equations 
and is described by a coupled system of random Schrodinger equa- 
tions driven by a Brownian field whose covariance is determined by 
the two-point statistics of the fluctuations of the random medium. 
For the reflected and transmitted fields the associated Wigner distri- 
butions and the autocorrelation functions are determined by a closed 
system of transport equations. The Wigner distribution is then used 
to describe the enhanced backscattering phenomenon for the reflected 
field. 

1. Introduction. The paraxial wave equation, in homogeneous or in ran- 
dom media, is a model used for many applications, for instance in commu- 
nication and imaging [19]. It has the form of an evolution equation that 
describes waves propagating along a privileged axis and it can be obtained 
by neglecting backscattering. Its simplicity, compared to the full three- 
dimensional wave equation, enables analysis of many important phenomena, 
such as laser beam propagation [9, 12, 16, 23], time reversal in random media 
[5, 11], underwater acoustics [24] or migration problems in geophysics [6]. 

The derivation of the paraxial model in homogeneous media is well un- 
derstood, and it can also be justified in heterogeneous media for small vari- 
ations of the wave speed [2, 4]. However, it is not clear whether the paraxial 
(parabolic) approximation is still valid in a scaling regime in which the 
medium fluctuations are rapid and can be approximated by a white-noise 
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term. The main motivations for studying the white-noise paraxial wave equa- 
tion are (i) it appears as a very natural model in many applications where 
the correlation length of the medium is relatively small, in particular much 
smaller than the propagation distance, (ii) it allows for the use of Ito's 
stochastic calculus, which in turn enables the closure of the hierarchy of 
moment equations [15] and thereby analysis of important wave propagation 
problems, such as the star scintillation due to atmospheric turbulence [25]. 

If the paraxial approximation and the white-noise approximation can be 
justified simultaneously, then the conjecture [21, 23] is that the limit equa- 
tion should take the form of the random Schrodinger equation studied in 
particular in [8]. The proof of the convergence of the solution of the wave 
equation in random media to the solution of the white-noise paraxial equa- 
tion was obtained in the case of stratified weakly fluctuating media in [1]. 
In our paper we consider the transmission and reflection of acoustic waves 
by a slab of medium whose parameters have three-dimensional random fluc- 
tuations and whose end is either transparent or a strong interface. This 
model is particularly interesting in the context of optical coherence tomog- 
raphy [26]. We analyze a wave propagation regime in which the paraxial 
and white-noise approximations are valid. In this regime we obtain a sys- 
tem of coupled Schrodinger equations driven by a Brownian field that fully 
determines the statistics of the transmitted and reflected waves. As a corol- 
lary we compute explicitly the two-point statistics of the transmitted and 
reflected waves. These results show that the often used "independent ap- 
proach" for the reflected wave (in which the statistics of the forward- and 
backward-propagating waves are assumed to be independent [26]) is valid if 
and only if the transverse correlation radius of the fluctuations of the random 
medium is smaller than the initial beam width. Finally, we use the coupled 
Schrodinger equations to give a rigorous account for the enhanced backscat- 
tering or weak localization phenomenon [27] and we compute explicitly the 
enhancement factor and the shape of the enhanced backscattering cone. 

2. The transmission and reflection operators. We consider linear acous- 
tic waves propagating in 1 -|- d spatial dimensions with random medium 
fluctuations. The governing equations are 



where p is the pressure field, u is the velocity field, p is the density of the 
medium, K is the bulk modulus of the medium and {z,x.) G M x M*^ are 
the space coordinates. The source is modelled by the forcing term F. We 
consider in this paper the situation in which a random slab occupying the 
section z S (0, L) is sandwiched in between two homogeneous half-spaces. 
The source, F, is located outside of the slab at z = zq, zq > L. We shall 





1 dp 
K{z,x.)dt 
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refer to waves propagating in the direction with a positive z component 
as right-propagating waves. The medium fluctuations in the random slab 
(0, L) vary rapidly in space while the "background" medium is constant. The 
medium is assumed to be matched at the right boundary z = L. We consider 
a possible mismatch at the boundary 2 = and denote the background 
medium parameters by po and Kq in the half-space z < and by pi and Ki 
in the half-space z > 0: 



where the random field z^/^(z,x) models the medium fluctuations, whose 
correlation length is Ik- The source has the form 



where is the unit vector pointing in the z-direction, zq> L \s the source 
position. We denote by uj^ the typical frequency of the source term fs and 
by Rq the diameter of its spatial support (which gives the initial beam 
width). The typical wavelength associated with the typical frequency ujq is 
Ao = 2-kci/uq, for ci = ^K\l p\ the background speed for z > 0, which is 
of the same order as the background speed cq = \/ Kq / pQ in the half-space 



We can now introduce the scaling regime that we consider in this paper: 

(1) We assume that the correlation length Ik of the medium is much 
smaller than the propagation distance L. We denote by the ratio between 
the correlation length and the typical propagation distance. 

(2) We assume that the transverse width of the source i?o and the cor- 
relation length of the medium Ik are of the same order. This means that 
we assume that the ratio Rq/L is of order e^. This scaling is motivated by 
the fact that, in this regime, there is a nontrivial interaction between the 
fluctuations of the medium and the beam. 

(3) We assume that the typical wavelength Aq is much smaller than the 
propagation distance L; more precisely, we assume that the ratio Xq/L is of 
order e^. This high-frequency scaling is motivated by the following consid- 
erations. The Rayleigh length for a beam with initial width Rq and central 
wavelength Aq is of the order of Rq/Xq in absence of random fluctuations 
(the Rayleigh length is the distance from beam waist where the beam area is 
doubled by diffraction). In order to get a Rayleigh length of the order of the 
propagation distance L, the ratio Xq/L must be of order since Rq/L ^ . 




if z < 0, 
ff zG (0,L) 
ff z > L, 




ff z < 0, 
ff zG (0,L) 
ff z > L, 



F(t,z,x) = fs{t,:x.)5{z - zo)e. 



■Z ) 



z<0. 
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Henceforth we shall assume nondimensionalized units chosen such that the 
background bulk modulus Ki and density pi in the half-space z > are 1, 
hence, the background speed ci = ^/K^Jpi and impedance Zi = ^/Kipi are 
also equal to 1. If we consider the propagation distance, L, as our reference 
distance of order 1 in this scaled regime, then the source has the form 

t 



(2.2) 



F(t,z,x) = / 



5{z - zqY 



where f{t,x) is the normalized source shape function (with time and spatial 
scales of variations of order 1), and the medium fluctuations have the form 

if z < 0, 

if zG (0,L), 

if z > L, 



K(z,x) 




p{z,x) 



if z < 0, 
if zG (0,L), 
if z > L, 

where the zero-mean, stationary random field has a correlation length of 
order 1 and standard deviation of order 1. We also assume that it satisfies 
strong mixing conditions in z. Here the amplitude of the fluctuations 
has been chosen so as to obtain an effective regime of order 1 when e goes 
to zero. That is, if the magnitude of the fluctuations is smaller than e^, 
then the wave would propagate as if the medium were homogeneous, while 
if the order of magnitude is larger, then the wave would not penetrate the 
slab. The scaling that we consider here corresponds to the physically most 
interesting situation. 

Since both the medium and the source have transverse spatial variations 
at the scale e^, it is convenient to rescale the transverse variable x/e^ — > x 
and to introduce the rescaled fields and p^: 



,z,e X . 



(2.3) u^(t,z,x) = u(t,z,e x), p%t, z,x) = p{t, 

The reader should keep in mind that thus, in the discussion below, when 
we refer to the transversal spatial parameter x it corresponds to e^x in the 
original coordinates. The rescaled fields satisfy in the region z G (— oo,0]: 



(2.4) po 



7U" 

'dt 



+ 



u 



0, 



1 dp^ 

where Vx stands for the gradient with respect to the transverse spatial 
variables x. In the random slab z G (0, L) the fields satisfy 



(2.5) 



~dt 



+ 



dz 



/ = 0, 



,x 



dp'^ 
'dt 



+ 



dz 



• u 
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Fig. 1. Boundary conditions for the modes in the presence of an interface at z — 0, a 
random slab (0,L) and a source at z = zq. 



and in the region z G [L^co) (in which the source is located): 



+ 



P 



~dt 



+ 



u 



0. 



In the Fourier domain (with respect to time), the first equations of these 
systems give the expressions of the velocity fields in terms of the pressure 
fields. 

We first consider the wave field in the homogeneous half-space z < 0, 
which allows us to introduce the standard paraxial wave equation in ho- 
mogeneous medium. The wave field satisfies (2.4) with the wave speed 
Co = v^o/pq. Following [18] (in a different scaling) we introduce the com- 
plex amplitudes Og and fog of the right- and left-propagating modes (see 
Figure 1): 



ao{k,z,x) = — 



bfi{k,z,x) = y 



J.kz/e'^ 



;/(l/(M,x) + if(M,x)).--/«** 

They are such that the pressure field in the region z < can be written as 

p'{t, z,x) = — f {a'oik, z, x)e*^"/"' + b'oik, z, ^y-ik-le^y-^'"-'k'l'-^ dk, 
2tt J 



and they satisfy 



5^(fc,.,x)e^^^/^Vp(A;,z,x)e-^^/^^ 
oz oz 



0. 



Using (2.4), we find that they also satisfy the coupled mode equations 

where Ax is the transverse Laplacian. In the limit e — > 0, the cross terms 
(proportional to 6=*=^*^^/^ ) average out to zero and we get the two uncoupled 
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paraxial wave equations 
(9ag _ 



dz 2k °' dz 2A; ° 
Taking into account the fact that there is no source in the half-space z <0, 
and therefore no right-going wave, we obtain the general expression of the 
wave in the left homogeneous half-space 



(2.6) p' 



-it,z,^) = ^ fmk,z,^)e-''''/'\-''°'''/''dk, z<0. 
2ir J 



Similarly, the wave fields in the homogeneous regions {L,zq) and {zo,oo) 
have the forms 



/(t,z,x) 



^ J ial{k,z,^)e'^'/'' +bl{k,z,^)e-''''/'')e-''''/'' dk, zG{L,zo), 
z,^y^^l^\-^^^l^^ dk, 



1 

2^ 



Z > Zq, 



respectively. Here we used the fact that there is no source and therefore no 
left-going wave in the region z > zq; see Figure 1. We can also use the jump 
conditions across the source position z = zq to obtain the relations 



(2.7) 
(2.8) 



bl{k,zo,x) 



2 



e—"'- /(A:,x), 



a|(/i;,zo,x) -al{k,ZQ,x) = \e 



1 -ifczo/e*^ 



/>,x). 



By solving the paraxial wave equation for 6f , we obtain the expression for 
the complex amplitude of the wave impinging on the random slab at z = L: 

(2.9) 6f(fc,L,x) = e^'=^«/^'6inc(fc,x), 

(2.10) 6inc(fc, x) = / f{K ,^)e^l-l'(^--o)/{2fc)+i«.x 

where the Fourier transforms are defined by 
fit,^)e''''dt, f{k,K)-- 



(2.11) /(A;,x 
he pressu 
p''(i,z,x 



/(fc,x)e-^'^-^(ix. 

The pressure field in the region z G (0, L) can be written as 
1 



2tt 



{d^{k,z,^)e'^^l^ +b%k,z,x.)e-''''/^ )e-iW^ dk, 



with the complex amplitudes and b^ of the right- and left-propagating 
modes given explicitly by 



a'^(/c, z,x) 
b^{k, z,x) 



1 1 dp'^ 

1 1 dp^ 



(t,z,x)^e^^*/"' dt 
e*^'*/^' dt 



-ikz/e'^ 
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daf ( ik ( z \ i ^ 



(2.12) 



(2.13) 



dz \2e Ve^' / 2k 

dz \2e W J 2k 

f ik f z \ i ^ 



This system is vaUd in z € (0, L) and it is complemented with the following 
boundary conditions at z = and z = L: 

(2.14) h'{k, z = L,x)= e"'''/'%inc(.k,x), 

(2.15) d'ik, z = 0, x) = TZob^ik, z = 0,x), 

where TZq = {Zq — 1)/{Zq + 1) is the reflection coefficient of the interface at 
z = and Zq = \JKqpq is the impedance of the left homogeneous half-space. 
These boundary conditions are obtained from the continuity relations of the 
fields p*" and 6;^ • at z = and z = L. The continuity relations also give 
the expressions for the complex amplitudes of the transmitted field in the 
region z < and of the reflected field of in the region z> L: 

(2.16) bl{k,z = 0,x)=Tob%k,z = 0,x), 

(2.17) af(/c,z = L,x) = a^(A;,z = L,x), 

1 /2 

where Tq = 2Zq /(I + Zq) is the transmission coefficient of the interface at 
z = 0. If there is no impedance contrast Zq = 1, then Tq = 1 and TZq = and 
the boundary condition (2.15) reads d^{k, z = 0,x) = 0. This is the radiation 
condition expressing the fact that there is no wave incoming from —oo. 

We now make use of an invariant imbedding step and introduce transmis- 
sion and reflection operators. First, we define the lateral Fourier modes 



a'^{k,z,K,)= / a^(A;,2;,x)e '^'^''^dx, 

(2.18) 

W{k,z,K)= b''{k,z,x)e-"^-''dx, 

and make the ansatz 

bl{k,0,K)= / T%k,z,K,K')b''{k,z,K')dK', 

d^{k,z,K) = / TV{k,z,K,K!)lf{k,z,K)dK. 
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Using the mode coupling equations (2.12)-(2.13) we find that the operators 
and IV^ satisfy 

dz 

(2.20) =e-2^'=^/^'£=(fc,z,K,K') 

^g2ifc^/.4 ffTie(^k,Z,K,Ki)C%k,Z,Ki,K2) 



X TV{k, Z, K2,k!) dKl dK2 
+ J C'^{k,Z,K,Ki)TZ'^{k,Z,Ki,K')dKi 

+ J TZ-^ {k, z, K, Ki)C'^ {k, z, Ki, k') dt^i. 



— f'{k,Z,K,K') 



(2.211 



J {k, z, K, Ki)C^ {k, z, Ki, k') dKl 

+ e^'^'/'' I j f%k, Z, K, Ki)C%k, Z, Ki,K2) 
X 'Rf[k, Z, K,2, k') dKl dK2, 



where we have defined 

(2.22) £^(A;,z,Ki,fv2) = - — \ki\'^5{ki - K2) + ^^^^^^^ t>(^^, - ^2 

with v{z,k) the partial Fourier transform (in x) of z^(z,x). This system is 
complemented with the initial conditions at 2 = 0, which are obtained from 
(2.15) and (2.16): 

n^k, z = Q,K, k') = no5{K - k'), f^k, z = 0,K, k') = %5{k - k'). 

The transmission and reflection operators evaluated at z = L carry all the 
relevant information about the random medium from the point of view of the 
transmitted and reflected waves, which are our main quantities of interest. 

Our objective in the next sections is to characterize the transmitted wave 
field 

Ptr('S)X) =p^(zo + e''s,^ = 0,x) 

(2.23) ^ 



(2vr) 



d+l 



T%k, L, K, K')binc{k, k') dK' e^^**-^-^^) dKdk, 
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and the reflected wave field 

Prof('S,x) =p'^{zo + L + e'^s,z = L,x.) 

n^k,L,K, K')binc {k, K.') dn e*('*-^-'='^) dn dk. 



(2.24) ^ 



(27r)^+i 

Note that the wave field "fronts" are observed on the time scale of the source 
and around their respective expected arrival times {zq for the transmitted 
wave, and L + zq for the refiected wave, which corresponds to the travel 
time to go from the source at z = zq to the interface at z = and back at 
the surface at z = L). 

3. The random Schrodinger model. 

3.1. Statement of the main result. We consider the transmitted and re- 
fiected fields pIj. and p^^j defined by (2.24)-(2.25) and use diffusion approx- 
imation theorems to identify a limit random Schrodinger model. The main 
result is the following one. 

Proposition 3.1. The processes (pjj.(s,x),p^j,f (s,x))gg]g ^g^d converge 
in distribution as e ^ in the space C^{R,L'^{M.'^,M?)) D L'^{R,L'^{M.'^,M.'^)) 
to the limit process (ptr(s,x),Prcf (s,x))sgiR^xGK'* 

(3.1) ptr{s,^) = ^ JJ t{k,L,^,^')kr.cik,^')d^' e-"'' dk, 

(3.2) p,ef(s,x) = i-yy7t(A;,L,x,x')6inc(A;,x')dx'e-^*^^dA:. 

Here C'^(M, L^(]R°',]R^)) is the space of continuous functions (in s) with 
values in L'^{R'^,R'^) and L"^ {R, L'^ {W^ ,R'^)) = L^{R x M'^,M2). The opera- 
tors 'T{k, z,x,x.') and 'R-{k, z,x,x.') are the solutions of the following Itd- 
Schrddinger diffusion models: 

(3.3) d'T{k,z,x,x) = — Ax'T(A;,2;,x,x')(iz + — T(A;,z,x,x') o dB(z,x), 

2k 2 

d'R,{k, z, X, x') = — (Ax + Ax')7?.(/c, z, x, x') dz 

(3.4) 

+ y7t(/c,z,x,x') o {dB{z,x) + dB{z,x)), 

with the initial conditions 

t{k, 0, X, X ) = To5{x - X ), Ti{k, 0, x, x ) = 7^o5(x - x ). 
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The symbol o stands for the Stratonovich stochastic integral in z, B{z,x.) is 
a real-valued Brownian field with covariance 

(3.5) E[B(zi,xi)B(z2,X2)] = min{zi,Z2}Co(xi - X2), 
and we have used the notation 

(3.6) C(z,x) =E[z^(z +z,x' + x)i/(z',x')], 



(3.7) 



Co(x) 



C{z, x) dz. 



The moments of the finite- dimensional distributions also converge, in the 
sense that 



E 



(3.8) 



n ^ ^jY"' n Prciisj , y 



£^0 



E 



J|ptr(Sj,Xj)"'^ ]^Pref(Si,Xi)' 

■j=l i=i 



for any G N, (sj)j=i,...,g G M", (sj)i=i,...,g S (x,),=i,...,g e M'^'?, 
(ij)j=i,...,gelR'^^ (m'i)j=i,...,<7GN'' and {mj)j=i,...,q . 

In [8] the existence and uniqueness have been estabhshed for the random 
process 



Vfc(z,x) = J T{k,z,x,x')(l){x')dx', 



for a test function cj) with unit L^(R'^)-norm, in the case Tq = 1. It is shown 
that the process Vk{z,x) is a continuous Markov diffusion process on the 
unit ball of L^(]R'^,C). The moment equations moreover satisfy a closed 
system at each order [15]. The analysis can be readily extended to the pair 
(Ptr;Pref) defined in terms of ('2", 7?.) and can be carried out jointly for all 
frequencies k in an interval bounded away from and infinity. We then get 
that the processes ptr and p^ef have constant L^(M x M'^)-norms, so that the 
conservation of energy relation holds: 

(3.9) JJ |ptr(s,x)p + |pi.ef(s,x)p(is(ix = yy |6inc(s,x)p(is(ix, 

where we have also used the identity TZq + Tq = I. 

The main steps of the proof of Proposition 3.1 given in the next section 
are 



(1) tightness and a priori estimates for (Ptr>Pref)' 
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(2) convergence of the finite-dimensional distributions of the process (pfj,, 
p^gj), using the convergence of specific moments of the reflection and trans- 
mission operators and 'K? . Heuristically, the reasons for the convergence 
of '7''^ and Tif are that the terms with the rapid phases exp(ibi2/i;z/e^) in 
(2.20)-(2.21) vanish in the limit e — > 0, and that the random "potential" 
in (2.22) can be replaced by a white noise. We then get formally the limit 
system (3.3)-(3.4). However, this holds true in a special weak sense only. In- 
deed, it is important to note that the reflection and transmission operators 

and 'K,^ themselves do not converge to T and TL solution of (3.3)- 
(3.4), but only certain moments (expectations of products of components 
with distinct frequencies k\ which are those needed to ensure the conver- 
gence of the transmitted and reflected fields. This approach is similar to the 
one used to prove the O'Doherty-Anstey theory in one-dimensional random 
media [7, 14] and to study the second-order statistics of the wave backscat- 
tered by a three-dimensional random medium [17, 18]. The limit moments 
are characterized by the system (A. 15) of coupled equations, which gives the 
practical way to compute all the moments of the reflected and transmitted 
wave fields. Some particular applications will be given in Sections 4 and 5. 

(3) use of Ito's lemma for Hilbert-space- valued processes [20], Theorem 
2.4, in order to check that the specific moments of the reflection and trans- 
mission operators 7" and 7?. given by (3.3)-(3.4) are solutions of the system 
(A. 15) of coupled equations. 

3.2. Proof of Proposition 3.1. This section is devoted to the proof of 
Proposition 3.1. We shall use a technique similar to the one presented in 
[14] in the case of randomly layered media. 

Step 1. A priori estimates. From (2.12)-(2.13) we can check that, for 
any /c, the integral 

j |a^(fc,z,x)|2- |6^(A:,z,x)|2(ix 

is conserved in z. Applying this conservation relation at z = and z = L, 
and taking into account the boundary conditions (2.14)-(2.15), we obtain 

j |a^(A:,L,x)|2(ix + (l-7^g) j |6-(A;, 0, x)!^ dx = j |6inc(A;,x)|2 dx. 

Using now (2.16)-(2.17) and the identity 7^g + T§ = 1, we obtain 

(3.10) j |a^(/c,L,x)|2(ix + j |6g(A;,0,x)|2dx = j |&inc(A;, x)!^ dx, 

which expresses the fact that the power of the incoming wave is fully re- 
covered by the transmitted and reflected waves. Integrating in k and using 
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Parseval's equality establishes the total energy conservation relation 
(3.11) JJ |p^ef(s,x)|2 (ix(is + JJ |pe^(s^x)p(ixds = JJ \binc{s,x)\'^ dxds. 
We first state a priori estimates for our quantities of interest. 

Lemma 3.1. There exists C > such that, uniformly in e and in so,si, 

J \pI^{sq,x)\'^ dx < C and 

J IPtr(siiX) -Ptr('50,x)p(ix < C\si - So\. 

The same estimate holds true for p^^^ . 



(3.12) 



Proof. Using Sobolev's embedding L°°(R) C i7^(M), there exists a con- 
stant Csob such that, for any x, 

C'sob 



sup |Ptr (S > X) I < Csob \\pIA-,^)\\h^ 



2-K 



'l + k^)\W{k,Q,x)\^dk, 



where we have also used Parseval's equality. Integrating in x and using the 
conservation equation (3.10) yields the first result of the lemma: 

Csob 



sup / \pI^{s,x)\ dx< / sup|ptr(s,x)| dx< 

s J J s 

By the Cauchy-Schwarz inequality, we have 



2lT 



'l + k^)\binc{k,x)\Uxdk. 



\pI,{si,x) -pI{so,x) 



Sl 



so 



< |si - Sol 



fSi rsi 

< ds 
J So J So 

2 



ds 



ds 



ds. 



The integral in x of the last term in the inequality can be bounded uni- 
formly as above. The reflected field can be analyzed in the same way, which 
completes the proof. □ 

Step 2. The moments of the finite- dimensional distribution of (]9jj.(s,x), 
p^g£(s,x)) converge to those of (ptr('S,x),prcf (sjx)). The general moment 
(3.8) of pfj.(s,x) can be expressed as the multiple integral 



E 



9 <? 

ij=i j=i 
1 



(27r)(^+^)('^+i) 
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X 



/■■■/ 



Y[ n '^'^'hj '^'^hj dkh,j 

h=lj=l 



^ n n '^^'hj di^hj dkh,j 

h=lj=l 

h,j 
h,j 



X E 



Lh,j h,j 



for N = X]ft=i "^/i M = Yl'h=i'^h- Therefore, the convergence of the 
general moment of the transmitted and reflected wave fields in the white- 
noise limit will follow from the convergence of the following specific moments 
E[P(L)] of the transmission and reflection operators, where 

N M 

(3.13) = \{f'{k.j,L,K.j,K'^)\{n'{kj,L,kj,k'^). 

We call these moments "specific" because we restrict our attention to the 
case in which the frequencies kj, kj are all distinct. 

We use diffusion approximation theorems to obtain equations for the mo- 
ments in the limit e — > 0. In the Appendix we show that 



limE[r(L)l =E 



N M 

IT '^(^j ■>L,Kj, Kj) 'TZ{kj,L, kj, kj 

Lj=l j=l 



when the right-hand side expectation is taken with respect to the following 
Ito-Schrodinger model for the transmission and reflection operators: 



d'T{k, z, K, k') 



l\2 



(3.14) 



+ 



2k 
ik 



T(k, z, K, k) dz ■ 



A:2Co(0) 



T{k^ z^ K, k) dz 



d'R-{k, z, K, k) 



2{2'kY 

i(|Kp + |Kf ) 

2k 

fe2Co(0) 



T{k, z, K, Ki) dB{z, Ki — k) dKi, 
'R,{k, z, K, k') dz 



7?.(A;, z, K, k') dz 
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(3.15) - J Co{Ki)'R.{k, z,K- Ki, k' - Ki) dni dz 

ik f 

+ T!.{k, z, K,i,K,') dB{z, n — Ki)) dKi, 

with the initial conditions 'T{k,0, k, k') =Tq5{k — k') and TZ{k,0, k, k') = 
TZq5{k — k'). Here we have used the notations (3.6)-(3.7) and the Brownian 
field B has the following operator-valued spatial covariance: 

(3.16) E[B(zi,fti)B(z2,K2)] =min{zi,Z2}(27r)'^Co(fti)5(Ki + K2), 
where 

/oo r 
/ C(z,x)e"''*-''(ix(iz. 
-oo Jm.<i 

The field B is the partial Fourier transform of the field B defined in the 
statement of the proposition. Consider next the reflection operator in the 
original spatial variables: 

(3.18) r(fc,z,x,xO = J j e'^''-''~'^'-''")r{k,z,K,K')dKdK', 

(3.19) 7t(A;,z,x,x') = J j e'^''-''^ ^ n{k, z, k, h,') dndn' . 

Then we find that this operator is weakly characterized by the following 
Ito-Schrodinger diffusion: 

d'T{k, z, X, x') = — Ax''7'(A;, z, x, x') dz ^ "T^k, z, x, x') dz 

+ yr(A;,z,x,x')rf5(z,x'), 

d'R,{k, z, x,x') = — (Ax + Ax') 7?.(/c, 2;,x, x') dz 
2k 

+ ^Tt{k,z,x,x'){dB{z,x) + dB{z,x')). 

In Stratonovich form this diffusion model becomes (3.3). This proves the 
last statement of the proposition (the convergence of the moments). 
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Step 3. Convergence of (pfr^Pref) (ptr,Prof) in C°(M, n 
L^(M, L^(M'^,R^)). Here is the space equipped with the weak topol- 
ogy. Lemma 3.1 shows that the process {plr^Prd) tight in C°(]R,L^(M'^,M^)). 
Moreover, the first estimate in the lemma shows that, for any L?^ 
functions (p^ip-, the random processes 

^lis) = / P?r(«>x)(/)(x)dx, 



y^{s) = j p^rf(s,x)V'(x) dx 



are uniformly bounded. Therefore, the finite-dimensional distributions are 
characterized by the moments of the form 



E 



where G N, nij^rhj G N, Sj,Sj G M, (l^j^ijjj G L^(M'^, 
can be written as multiple integrals: 



). These moments 



E 



■i=i 



1 



(27r)(^+^^)(^+i) 



q mh 



Yl Y[dK'f^jdKh,jdkh,j X n n dk'f^jdKhjdkhj 

h=lj=l h=lj=l 



X E 



h,j 



for N = J2l^imh and M = J2l^irhh- Note that only specific moments of 
quantities of the form (3.13) appear (i.e., moments of products of the trans- 
mission and refiection operators at distinct k). The convergence of these spe- 
cific moments therefore implies the convergence of the finite-dimensional dis- 
tributions, hence the weak convergence in 

C°(M,L2 (M'^,M2))_ Furthermore, 
the estimate (3.11) shows that the processes are tight in L^(M, L^(M'^,]R2)) 
(the unit ball is compact in the weak topology). This proves the weak con- 
vergence in (M,l2 (R'^,m2)). 
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Step 4. Convergence of {pl„p%f) to (ptr,Pref) in {R, {W^ ,M.'^)) n L"^ {R, 
L^(M'^,M^)). Here is the space equipped with the strong topology. 
Since the convergence has been proved in the weak topology, it is sufficient 
to show that the L^-norm is preserved. On the one hand, the L^-norms of the 
processes (pfi-jPref) deterministic, independent of e, and given by (3.11). 
On the other hand, the limit process {ptr,Prc{) has constant L^-norm given 
by (3.9). This proves the first statement of the proposition and completes 
its proof. 



C(z,x) 



4. The Wigner distributions. We first introduce the dimensionless au- 
tocorrelation function C of the fiuctuations of the random medium 

z X ' 

where a is the standard deviation of the fluctuations of the random medium 
and Iz (resp. Ix) is the longitudinal (resp. transverse) correlation radius of 
the medium. With this representation we have 

Co(x) = a\Co (^f^ , ao(u) = a%ltCo{nlx). 

We assume next that the power spectral density Co(u) decays fast enough 
so that /|upCo(u)(iu is finite. This means that the autocorrelation func- 
tion Co(x) is at least twice differentiable at x = 0, which corresponds to a 
smooth random medium. For simplicity, we assume also that the random 
fluctuations are isotropic in the transverse directions, in the sense that the 
autocorrelation function Co(x) depends only on |x|. 



4.1. The Wigner distribution of the transmitted wave. We now consider 
two frequencies ki and k2 in a frequency band centered at k and we define 
the two-frequency Wigner distribution of the transmission operator by 



(4.1) 



-«(q-y+q'-y') 



X E 



#rx'+^ 



Vk 



n^2 



y_ 

2 



dydy'. 



Using the stochastic equation (3.3) and Ito's formula, we find that the 
Wigner distribution satisfies the closed system 
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+ 4(2^ y C'o(u)iy,,,, (^z, X, X , q, q - - + ^ j u 

starting from 

W^J^,,, (z = 0, X, x', q, q') = To^{A7T^k,k2/kY/^5{^ - ^')5{q + q'). 

It is possible to solve this system and to find an integral representation of the 
two- frequency Wigner distribution by using the approach of [9]. However, we 
aim at focusing on spatial aspects in the next sections, and we shall simplify 
the algebra by assuming that the bandwidth B of the incoming wave with 
carrier wavenumber ko is small. To describe this regime it is convenient to 
introduce 

where tq is the initial beam width, (3 describes the intensity of random scat- 
tering, while a and ao represent the intensities of diffraction on respectively 
the scales of the medium variations and the input beam. Note that these pa- 
rameters correspond to inverse Fresnel numbers (up to a factor 27r) relative 
to respectively the lateral medium correlation scale and the aperture; below 
we shall consider explicitly the case with small medium Fresnel number. We 
assume that the bandwidth B of the incoming wave is small in the sense 
that 

(4.3) B<.Bc, Bc:=komm{l,a^^,aQ^,(3-^). 

In this regime we can approximate the two-frequency Wigner distribution 
by its behavior at the carrier frequency. That is, if ki, /c2 lie in the spectrum 
of the incoming wave, the two-frequency Wigner distribution VF^ can 

be approximated by the simplified Wigner distribution W'^ that depends 
only on the carrier wavenumber ko and not on the lag k\ — k2- This Wigner 
distribution can be written in the form 

T^^(z,x,x',q,q')=To2(27r)'^W^f^,-,-,q/,,q^ 
where W"^ satisfies 

(4.4) ^ + • Vx'W^ = ^ / Co(u)[W^(q' - u) - dn, 

starting from yV"^(C = 0,x,x',q, q') = (^(x — x')(^(q + q'). We remark that 
these are in fact the classical equations of radiative transport for angularly 
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resolved wave energy density [22]. In [22] the radiative transfer equations 
are written in the standard time-dependent form, while (4.4) is written in 
a time-harmonic form in a reference frame moving with the background 
velocity along the z-axis corresponding to the form set forth in [10] . 

In this context we have that /3Co(0) is the total scattering cross-section, 
/3Co(-) the differential scattering cross-section describing coupling of modes 
depending on their relative propagation directions and (a/a;/ro)q' a trans- 
port or velocity vector. It is clear from this equation that diffractive effects 
(characterized by the term q' • Vx') are of order 1 if alx/r^ ~ 1 or equiva- 
lently korolx ~ L. In terms of Fresnel numbers, diffractive effects are of order 
1 if 

_ L 

OLe = 2 ~ 1' 

where we have defined the effective aperture by 



In Section 5 we shall see that a physically important and mathematically 
interesting regime corresponds to 

ao Oe ~ 1 ^ a, 

so that from the point of view of the "medium Fresnel" number we are in 
a Fraunhofer diffraction scaling, while from the point of view of the "source 
Fresnel" number we are in a Fresnel diffraction scaling and finally from the 
point of view of the effective Fresnel number, l/(27rae), we are in a general 
or scalar diffraction theory setup. 

By taking a Fourier transform in q' and x', we obtain a transport equation 
that can be integrated and we find the following integral representation for 



W^(C,x,x',q,q' 
(47r2a 



7-2 

(4_5) _ / / g-j(q'+q)-»7i-«(»'o(x'-x)/(QZa;)+qC)-T72 



This expression will be used in Section 5.1 to compute and discuss the two- 
point statistics of the transmitted field. 

4.2. The Wigner distribution for the reflected wave. We define the two- 
frequency Wigner distribution of the reflection operator by 

^fc1,fc2(^>x,x',q,q') 
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(4.6) 



-j(q-y+q'-y') 



X E 



TZi ki,z, 



-I 



Vk 



xlZiko.z 



Vk 



y 

X 

2 



Vk 



2 



dydy'. 



If the bandwidth of the incoming wave satisfies (4.3) and if ki,k2 He in the 
spectrum of the wave, then we find by using (3.4) that the two-frequency 
Wigner distribution W^^ can be approximated by the simphfied Wigner 
distribution that depends only on the carrier wavenumber /cq a-nd not 
on the lag ki — k2- This Wigner distribution satisfies the closed system 

dz 



ko ko 

Kq 



4(27r)'^ 



Co{u) 



^-^(2;, X, x', q - u, q') + 1^-^(2:, X, x', q, q' - u) 



+ 2W^« ( z, X, x', q - ^u, q' - ) cos (u • (x - x')) 



2W 



R 



z,x,x ,q- ^u,q' + ^u) cos(u- (x-x)) 



2VF^(z,x,x',q,q') 



starting from W^{z = 0,x,x',q,q') = 7?.Q(27r)'^(5(x - x')(5(q + q'). We now 
cast the Wigner distribution in a suitable dimensionless form. We consider 
the following Fourier transform of the Wigner distribution W^: 



^(z,x,x',q,qO = -^/ y«(z,^,q-q',s)e-('''--)cis, 



w 



which we introduce because the stationary maps that we will identify in 
Lemma 4.1, in the asymptotic regime a ^ 00, have simple representations 
in this new frame. Note also that this ansatz incorporates the fact that 
does not depend on x + x', only on x — x', q and q', which follows from 
the stationarity of the random medium. The Fourier-transformed operator 
V^{z,q^,r,s) has the form 

V%, q, r, s) = 7^2(^/,)'^e^--/'=o rl,, sl^ , 
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where is the solution of the dimensionless system 



X / Co(u) 



+ V^(C,q-^u,r + u,s)e™^ 
(4.7) +V^(C,q- -u,r,s-u)e-™^ 



2 

+ (C, q - ^u, r, s + u) e*-"^ - 2V^(C, q, r, s) 



(C' q - ^U, r - U, S + gia[(r-s).u-|u|2]C 

_ V« ("c, q - iu, r - u, s - e--[(r+«)-+l"l^l?' 



du. 



starting from V'^(C,q, r, s) = 5(q). The parameters a and (3 are given by 
(4.2). 

We want now to analyze the regime in which the transverse correlation 
length Ix of the medium is smaller than the beam width rg. More exactly, 
we assume from now on in this section that 

(1) Ix, which means that the transverse correlation length of the 
medium is small, 

(2) korolx ~ -L, which means that diffractive effects are of order 1. 

These two conditions are equivalent to ao ctg ~ 1 <C a. Note that in the 
previous section we established the fact that diffraction plays a role for the 
transmitted wave for a propagation distance L of the order of korol^, which 
is smaller than the usual Rayleigh length kor^. This is a well-known result 
[19], and we shall deduce it for the reflected field in the analytic framework 
that we have set forth. 

The rapid transverse variations regime is particularly interesting to study 
because has a multiscale behavior. In (4.7) this regime gives rise to 
rapid phases. The following proposition describes the asymptotic behavior 
of as a—>-oo. The presence of singular layers at r = and at s = re- 
quires particular attention and is responsible, for instance, for the enhanced 
backscattering phenomenon studied in Section 5.3. The situation with a 
large corresponds to a strong diffraction situation, at the scale of the lateral 
medium fluctuations. In general [part (1) in Lemma 4.1] the intensity of the 
reflection operator decays exponentially according to the parameter /3Co(0) 
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corresponding to the total scattering cross section. This decay foUows from 
a partial loss of coherence by random forward scattering. However, as artic- 
ulated in parts (2) and (3) of the lemma below, the coupling of wave modes 
depends on the full medium autocorrelation function if we look at nearby 
specular reflection or small spatial offset frequencies. This coupling will be 
important when we analyze enhanced backscattering in Section 5.3. 

Lemma 4.1. (1) For any r^O, s^O: 
(4.8) V«(C, q, r, s) 5(q)e-2/5'^o 



(2) For any s^O we haveV^{C,q,^,s)°'-^^ V^{C,q) where V^{C,q) is 
solution of 



dV^ 213 



and is given explicitly by 



( C, q - ) cos(r • uC) - Vr^(C, q) 



1 



Rt 



(4.10) Vr^(C,q) = /e-''i-e^/o''^»("/2+<)+Co{u/2-rC')-2Co(o)dC' 
(2-kY J 

Similarly, for any r we have V^(C, q^ r, ^) '^H^ V^(C, q) ■ 
(3) For any r and s we have 



(4.11) C,q,-,- 

a a 



r s 



V«(C,q)+K^(C,q)-5(q)e 



~2/3Co(0)C 



Proof. In case (1), the rapid phases cancel the contributions of all but 
the term V"^ (C,q, r,s) in (4.7), and we get 



'^TTVd /Co(u)V^du = -2/3Co(0)V«, 



which gives (4.8). In case (2), we obtain in the limit a — > oo the simplified 
system 



R 



Co(u) 



Vr^(C,q-^u,s-u)e 



1 



+ Vr"" ( C, q - 2 s + u ) e^-"C - 2V«(C, q, s) 



du. 



We then Fourier transform this equation in q and s, and we obtain that the 
solution does not depend on s, that it satisfies (4.9), and that it is given by 
(4.10). 
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In case (3), we obtain the simplified system for V/^s(C, q) = limo^oo V'^(C, q, 



L. 

a'' a' 



dVr% 2(3 



dC (27r) 



Co(u) 



V,^(C,q-^u)cos(s-uC) 



1 



1 



+ ( C, q - 2" ) cos(r • uC) - Vr^,(C, q) 



Using (4.9) satisfied by and V^, we get 



which yields (4.11). □ 



dC 



du. 



5. Two-point statistics of the transmitted and reflected fields. 



5.1. The transmitted field. The results of Section 4.1 allow us to compute 
the two-point statistics of the transmitted field. We assume that 

(a) the pulse has carrier frequency /cq and it is narrowband in the sense 
that it satisfies (4.3), 

(b) the input beam spatial profile is Gaussian with radius ro, 



(5.1) 



6inc(t,x) = fo{t)e 



-ikot 



exp 



(suppressing the complex conjugate part here and below), 

(c) the transverse correlation radius Ix of the random fluctuations of the 
medium is much smaller than ro and korolx ~ L. As we have discussed after 
(4.4), this last condition ensures that diffractive effects are of order 1. 

Under (a) and (b), we find that the autocorrelation function of the trans- 
mitted field defined by 



^tr(s,*,x,y) = liniE 



has the form 



Ptv[ ^+ 2'^"^ 2 



t\^ 



Ptr 1 s 



t y 

— , X 

2' 2 



At,{s, t, X, y) = %ffo {s + -]fo[s--] e-'' 



ikot f ^0 \ 



2 \ d/2 



(5.2) 



-|r,L/fco+y|V(2r-^)-'-ol^lV8p-ir,-x 



X e 



kl/A J^Coivz/ko+y)-Co{0)dz 



drj. 
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Under (a)-(c), we obtain 

^tr(s 

(5.3) 



A,{s, t, X, y) = ToVo {^ + 1)To(s-1) e"^^«* (^) 



rg|r7|V8-|y|V(2rg)g-ir,.Xgfcg/4/^ Co(riz/ko+y)-CoiO) dz _ 

If, moreover, random scattering is strong, in the sense that /3 ^ 1, or equiv- 
alently A;qCo(0)L ^ 1, and if the autocorrelation function of the random 
fluctuations of the medium is twice differentiable at zero: 

Co(x)=Co(0)-:^|x|2 + o(|x|2), 

(5.4) 

I^ = -^ACo(0) = -^ACo(0), 
then we obtain that the autocorrelation function has the Gaussian shape 

i\ — f. t\^~ikot( 



Atr(s,t,x,y) fo(s + - ]fo(s-- ]e 

(5.5) 



, 2|x|2 |yP , . x-y 

X exp — — ^ - - — —J + 1- 



The beam radius rx{L), the correlation radius Pt{L) and the parameter 
Xt{L) are characterized by 



(5.6) ^^(L)=roJl + ^, 



3r^ 



(5.7) pt{L) = ro 

(5.8) xt{L) 



l + L»L3/(3r, 



2^ 



1 + klrlDL/A + 
rriL) 



where we have taken into account that kor^ S> -L in the considered regime. 
Note in particular that the beam width increases at the anomalous rate 
L3/2 (which was first obtained in the physical literature in [13] and con- 
firmed mathematically in [11]). Furthermore, the lateral correlation radius 
decays to zero, which means that the beam becomes partially coherent. We 
finally remark that these results hold true in the case with a smooth random 
medium [with Cq twice differentiable at as in (5.4)]; the situation with a 
rough random medium will be addressed elsewhere. 
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5.2. The reflected field. Under (a) and (b), the limit autocorrelation 
function of the reflected field defined by 



Ar.ei{s,t,x,y) = limE 

has the form 

-4ref(s,i,x,y) 



Prcf(s+2'^+ 2 j^r^fr ~ 2'^ 2 



(5.9) 

dri^ dri2 (iT73e*^^3-^2/'=o-^o{l^3p+l^2-2r7iP)/8 



where the function is the solution of the system (4.7). Under (a)-(c) the 
asymptotic behavior of the function is determined by Lemma 4.1 and 
we find that the autocorrelation is given by 

A,,{s, t, X, y) = nlfo + - e-'^o* (^^y 

(5.10) X y'e-^ol^lV8-|yP/(2r2)g-i^.x 

X e'^o/'^ J^'' Co(vz/ko+y)-CoiO)dz 

This is exactly the form (5.3) of the autocorrelation function of the trans- 
mitted wave, upon the substitution 2L for L. This shows that we would have 
obtained the same result if we had assumed that the backward propagation 
was independent of the forward propagation. The independent approach is 
valid in the regime in which the beam width is much larger than the trans- 
verse correlation radius of the fluctuations of the random medium, but it is 
not valid in other regimes, as can be seen by comparing the full expressions 
(5.2) and (5.9). 



5.3. Enhanced backscattering. The comparison of the autocorrelation func- 
tion of the reflected wave and that of the transmitted wave for the propaga- 
tion distance 2L shows that, in the regime a ^ 1, there is no coherent effect 
building up between the forward and backward propagations. However, cor- 
rective terms show that there are some residual effects. In particular, we 
would like to show that the reflected intensity exhibits a singular picture in 
a very narrow cone, of angular width of order , around the backscat- 
tered direction. This phenomenon, called enhanced backscattering or weak 
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localization, is widely discussed in the physical literature [3, 27]. The phys- 
ical observation is that, for an incoming quasi-monochromatic quasi-plane 
wave, the mean reflected power has a local maximum in the backscattered 
direction, which is twice as large as the mean reflected power in the other 
directions. 

In this section, we assume that the incoming wave has the form 

6inc(i,x)=/(t)e-^'=o*ffinc(x), 

that it is narrowband in the sense that it satisfies (4.3), and that it is nearly 
a plane wave, in the sense that ^inc(«^) is concentrated at some Kinc (assumed 
to be different from the normal incident vector 0). By "concentrated" we 
mean that the angular width of the incoming beam is smaller than a~^. The 
reflected signal in the direction kq is 

= ^J n%k, L, KO, K%^,ik, K')e-'''' dk. 

The moment of the square modulus of ^^^^(s, fto) only involves specific mo- 
ments of quantities of the form (3.13) (with distinct k). Therefore this mo- 
ment converges to the one of the limit process Pj-c{{s,kq) defined as the 
Fourier transform in x of Prcfis,x) given by (3.2). This means that the 
mean reflected intensity in the direction kq converges to 

EmAs,>^o)\']'-^nl\f{s)\'i^{Ko), 

Using the fact that ginc('^) is concentrated at ftmc, we get 

(5.11) I^iKo) = PV^ (l, '^^^^l., {ko + /.inc)/., O) , 

where P = 2'% J \gir,c{t^'i)\^ dK[. This formula gives the mean reflected in- 
tensity in the direction ko and is valid for arbitrary values of a and (3. Let 
us consider the regime a S> 1. The mean reflected intensity far enough from 
the backscattered direction — ftinc is of the form 

/«(ko) = PVo^(i,^^^^^/x) for \Ko + K,^,\l,»a~\ 

where we have used the second point of Lemma 4.1. In a narrow angular cone 
around the backscattered direction —Kinc the reflected intensity is locally 
larger: 

/^(-/^inc + a- V) = P[V^(1, -Kinc^x) + V^;^ (1, -Ki^J,)], 
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where we have used the third point of Lemma 4.1. If we assume, additionaUy, 
that 1, then we have 

where V is the dimensionless version of D given by (5.4): D = a'^l^l^'^'D. This 
formula gives the width of the diffusion cone around the specular direction 



(5.12) AKspec = ^^ = ^^ f^^ = ./DLkQ 



On the top of this broad cone, we have a narrow cone of relative maximum 
equal to 2 centered along the backscattered direction —i^inc'- 

+ a- V) = p(vrP/3)-<^/2e-l-incP^S/{2'/3) [i + 
This shows that the width of the enhanced backscattering cone is 



(5.13) Akebc 



Note that the angular width A^ebc = '^'^ebc / of the cone is proportional 
to the wavelength, as predicted by physical arguments based on diagram- 
matic expansions [27]. 

APPENDIX: DERIVATION OF THE LIMIT MOMENT EQUATIONS 

The purpose of this appendix is to compute the limit of the expectation 
of (3.13) as e — > 0, for distinct frequencies kj, kj. Using (2.20) and (2.21) we 
find 

— (z) = 5: n r^ki,Z,Ki,K[) 

j = l l=ljLj 

M 

xWTV{kuz,kuk'i) 

1=1 

i^J I' [kj , Z, Kj , [kj , Z, Ka, Kj) dKa 

-\- e ' ^ ^ J (^jf) ^) '^jri f^a)-^ {kj , Z, Ka, Kf)) 
X 'R.^{kj,Z, Kf), Kj) dKa dKf, 

M N 

, Z, Kl, Ki) 

j = ll = l 



X 



COUPLED PARAXIAL WAVE EQUATIONS IN RANDOM MEDIA 27 

(A.14) 



M ^ 



X 



7?. [kj , Z, Kj , Ka)^ [kj , Z, Ka, H-b) 
X TZ^{kj,Z, Ki), k'j) dKa dKh 
-\- J (kj , Z, kj , {kj , Z, Ka, kj) dKa 

+ ^7?. (kj , Z^kj , Ka)C (kj , Z, kj) dHa 

We next apply the diffusion approximation to get limit equations for the mo- 
ments; see [14] for background material on and related applications of the 
diffusion approximation. Observe that the random coefficients are rapidly 
fluctuating in view of (2.22). Those coefficients that are of order are 
centered and fluctuate on the scale e^; moreover they are assumed to be 
rapidly mixing, giving a white-noise scaling situation. Moreover, the rapid 
phase terms exp(ib2iA;2;/e^) lead to some cancellations between interacting 
terms. Here, the fact that the frequencies are distinct plays a key role. As 
a consequence, by applying diffusion approximation results, we obtain the 
equations for the moments E[P] in the limit e — > 0: 

7(z) = limE[r(z)]. 

We obtain from (A.14) that / solves a system of integro-differential equa- 
tions 

N M 



\j=l j=l / 



< kjkiI{K'j -k,k'i + k) 
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M 



(A.15) 



+ 2 ^ kjI{Kj — K, k'j — k) 

M 

+ ^^kjki{I{kj -K,k'i- k) 

+ I{ki - K, k'j - k) 
+ I{kj - K,ki + k) 
+ I{k'j - K,k'i + 

AT M 

+ 2^^kjki{I{K'j - K,ki -k) 
j=i 1=1 



+ I {k'j — K,k'i + k)) > dK, 



where we only write the shifted arguments for /. The initial conditions are 



I{k, k, K, k, k', k', z = 0)= riiLi Hi^j - '^'j)T^o Y{f=i 5{i^3 - f^'j)- Using in 



particular the relation 



E 



JO 



Aa(Sa, Ka)Xb{Sb, Kfe) dB{Sa, K^) dB{sb, Kfo) dl^a dn^ 

min{2a,Zi,) 



E [\a{s, K)Xb{s, -k)] (27r)'*Co(K) ds dK, 



we can then verify that 
I(z)=E 



' N ^ M 

■j=l 3=1 



where the right-hand side expectation is taken with respect to the Ito- 
Schrodinger model for the transmission and reflection operators in (3.14)- 
(3.15). 
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